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ABSTRACT 

The total extragalactic background radiation can be an important test of the global 
star formation history (SFH). Using direct observational estimates of the SFH, along 
with standard assumptions about the initial mass function (IMF), we calculate the 
total extragalactic background radiation and the observed stellar density today. We 
show that plausible SFHs allow a significant range in each quantity, but that their ratio 
is very tightly constrained. Current estimates of the stellar mass and extragalactic 
background are difficult to reconcile, as long as the IMF is fixed to the Salpeter slope 
above 1 Mq. The joint confidence interval of these two quantities only agrees with that 
determined from the allowed range of SFH fits at the 3a level, and for our best-fit values 
the discrepancy is about a factor of two. Alternative energy sources that contribute to 
the background, such as active galactic nuclei (AGN), Population III stars, or decaying 
particles, appear unlikely to resolve the discrepancy. However, changes to the IMF 
allow plausible solutions to the background problem. The simplest is an average IMF 
with an increased contribution from stars around 1.5-4 Mq. A "paunchy" IMF of 
this sort could emerge as a global average if low mass star formation is suppressed in 
galaxies experiencing rapid starbursts. Such an IMF is consistent with observations 
of star-forming regions, and would help to reconcile the fossil record of star formation 
with the directly observed SFH. 

Key words: galaxies: stellar content - cosmology: diffuse radiation 



1 INTRODUCTION 

In the last decade, observations have built up a picture of 
when and where stars form in the universe. From inventories 
of the stellar co ntent of the local universe llCole et al.ll200ll : 
iBell et al.ll2003ll and sur veys of star format i on over its en- 
tiro cosmic history (e.g. iMadau et al.lll99(i lOabasch et alJ 
12004), it appears that star formation was much more rapid 
in the past, reaching a gentle peak at 2; ~ 1-3 and falling off 
towards higher redshifts. These general trends are in accord 
with the A-dominated cold dark matter (LCDM) model of 
structure formation that has become widely accepted in re- 
cent years. This model successfully explains the spectrum of 
perturbations in the microwave background, the clustering 
properties of local galaxies, the optical depth and fluctuation 
spectrum of the Lyg forest ga s, and the flux from distant su- 
pernovae iSoergel et al ."2006). However, galaxy formation is 
less well understood than these other phenomena. Theoret- 
ical models of galaxy formation can be i mplemented either 
by u s ing semi-an alytic calculat ions (e.g. |KauffaiMLii £t_al| 
I1999I: ICole et aljEoOO.: .SomcrviUe. Primack. fc FabeJl200ll) 



or by direct numerical simulations (e.g. iKatz et al.l 11993 : 
[Sprinacl & HcruQuist 2003); each method has its own well- 
known advantages and drawbacks. Difficult issues that af- 
fect the results in both methods include galaxy supernova 
and A GN feedback dSpringel et al.ll2005l:lDekel fc BirnboimI 
I2OO6), the history and physical influence of reionisation 
jGnedin 2000 ; Bullock et al. 2000) , the influence of "cold 
mode" accretion JKeres et al.ll2005fl. and th e properties of 
zero-metallicity stars IITumlinson et alJl2004r) . Different as- 
sumptions about these issues can produce quite different 
histories of star formation, so figuring out the total amount 
of star formation is a vital clue to the physics. 

The observational measures of star formation are suffi- 
ciently imprecise that cross-checks using different measures 
are important. The star formation history (SFH) is mea- 
sured by the direct emission from young stars or repro- 
cessed radiation from dusty star-forming regions. This can 
be tested against the evolution of the stellar mass, using light 
from old stars (preferably in the near-infrared) as a proxy 
for mass; this has been carried out by many authors (e.g., 
ICole et all I2001I: iDickinson etal\ 120031: iGwvn fc Hartwickl 
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I2OO5I: iHopkins fc Beaco rn 2006). The integrated star for- 
mation rate tends to exceed the observed stellar mass, es- 
pecially at high redshift. Given the large statistical and 
calibration errors in both quantities, however, the evi- 
dence for a discrepancy is not ve ry strong. The observed 
metal density can be used as well (iMadau fc Pozzettill2000l : 
iHopkins fc Beacomll200 m. though the results from this ap- 
proach are not yet precise enou gh to indicate more tha n a 
very approximate agreement. IHopkins fc BeacomI (|2C)09^ in- 
troduced the use of neutrinos from Type II supernovae to 
place an upper limit on the normalisation of the star for- 
mation history. Uncertainties in the effective neutrino tem- 
perature and stellar mass thresholds of Type II SNe cause 
systematic uncertainties in this limit; with present data, it 
is probably best to view the limit as constraining the super- 
nova physics given the estimated SFH. 

The extragalactic background light constitutes another 
vital test of the star formation history, since it records the 
total stellar emission over all time, weighted by the scale fac- 
tor at the time of emission. The observed background light 
contains two humps, one in the far infrared (FIR) and one in 
the ultraviolet (UV) to near-infrared (NIR) range. The FIR 
hump primarily owes to thermal dust emission from heavily 
obscured, rapidly star-forming galaxies. There is much de- 
bate over the optical/NIR portion, in which the sum of the 
observed galaxies amounts to < 50% of the total detected 
backgr ound (Wright 2001; Bern stein. Freedman. fc Madord 
120028^ . The observed spectrum in this region might contain 
a hint of Lyman- a emis sion from low-metallicit y, high-mass 
Pop III stars at z > 6 JMatsumoto et alj|2005r . but obser- 
vational uncertainties and theoretical difficulties make this 
idea very controversial. It is certainly possible that most 
of the excess optical/NIR background light comes from the 
outer reg ions of normal galaxies or stars in faint undetected 
galaxies JBernstein. Freedman. fc Madorell2002lJ) . 

Many authors have attempted to tie the background 
radiation in different wavebands to the history of specific 
subsets of galaxies, e.g. FIR emission from dusty starbursts 
or the NIR emissi o n from pre-reionisation epoch galaxies. 
iMadau fc Pozzettil i200Ct ) were the first to consider the to- 
tal amount of the background radiation as a record of the 
cumulative stellar emission from all types of galaxies. This 
exercise is valuable because reprocessing of stellar radiation 
by dust and atomic continuum opacity, which can be very 
difficult to model correctly, will in general not significantly 
change the total amount of emitted energy. The stellar back- 
ground energy is dominated by the stars that have already 
evolved off the main sequence, down as low as 0.9 Mq, so 
it samples a broader range of stellar masses than the di- 
rect indicators_of_the SFH. Madau fc Pozzett i (2000; see 
also lMadau et al.l200lL|Pozzetti fc Madaul26oil) found good 
agreement between the background radiation and the total 
amount of stars in the universe, concluding that there was 
little room for alternative sources of energy such as an early 
generation of stars with a top-heavy IMF. Since the work 
of Madau fc Pozzetti, additional measurements of the back- 
ground have been made in several different bands, and our 
understanding of the star formation history and local stellar 
density has improved. It thus seems worthwhile to test the 
background radiation as an indicator of the star formation 
history once more. 

In f|^ we discuss how we calibrate various indicators of 



the stellar mass and the star formation rate, and discuss the 
range of masses contributing to each. In f|3 we discuss the 
observed cosmological background radiation, star formation 
history, and stellar mass density in the universe and we con- 
struct an ensemble of fits to the SFH to test the sensitivity of 
other observables to it. In [21 we compare the results of these 
various indicators of star formation to one another. Using 
the ensemble of star formation histories, we conclude that 
it is difficult to reconcile the background radiation with the 
stellar mass density observed today, using our default IMF 
assumptions. i|K| discusses several possible solutions to this 
problem, focusing on the role of the stellar IMF. We show 
that the observed ratio of the total background light to the 
present-day NIR luminosity density is a powerful test of the 
IMF. Despite the large observational uncertainties, we can 
significantly restrict the allowed set of IMFs. i|S| summarises 
our conclusions. 

Since the IMF figures in every stage of our analysis and 
emerges as a focal point of our discussion, it is worth making 
a few introductory comments about its role. Indicators of the 
instantaneous star formation rate are dominated by massive 
stars; in our standard calculations below, for example, half of 
the UV luminosity comes from stars with mass M > 15 Mq. 
The contribution of a stellar population to the extragalactic 
background depends on the integrated bolometric luminos- 
ity over its lifetime. This includes a significant contribution 
from less massive stars, and in our standard calculations half 
of the background comes from stars with mass M > SM© . 
The present if-band light density comes largely from sub- 
giants and giants with masses just above the main-sequence 
turnoff, typically M ~ IM© for an old stellar population. 
The expected ratio of these three measures of the cosmic 
star formation history, therefore, depends on the shape of 
the IMF above IMq . Much of the mass in the IMF resides 
in low mass stars, and there has been a great deal of obser- 
vational investigation of the IMF shape below IMq. How- 
ever, low mass stars emit very little light, so IMF changes 
in this regime have almost no impact on the relative values 
of these three probes of the cosmic star formation history; 
they change the expected/inferred stellar mass density for 
the three measures by the same factor. Our conclusions in 
this paper are, therefore, almost entirely insensitive to the 
adopted shape of the IMF below ~ O.SM©. 

Throughout this paper, we assume a flat "737" 
LCDM cosmology of Q.m = 0.3, 0.^ = 0.7, h = 
iyo/(100kms"^Mpc"^) = 0.7, and 0,^ = 0.02 ft"^ con- 
sistent with values derived from the WMAP experiment 
JSpergel et al.ll2006H . 



2 CALIBRATIONS OF STELLAR ENERGY 
OUTPUT 

In the central parts of some galaxies, the mass in stars can be 
estimated directly from its dynamical effects, but on larger 
scales the dynamics are dominated by dark matter. Hence, 
at present, all direct estimates of the global production of 
stars are based on the light they produce. The calibration of 
the mass-to- light ratio is sensitive to the age and metallicity 
of the population, but these can be constrained using the 
observed colours and/or spectral lines. The other major un- 
certainty is the stellar initial mass function (IMF). The most 
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Table 1. Initial Mass Functions 



Form" 


K 


X 


Ml 


Mu 


Salpctcr 


0.172 


1.35 


0.1 


100 


diet Salpeter 


0.218 


1.35 


0.188 


100 


Kennicutt 


0.328 


0.4 


0.1 


1 




0.328 


1.5 


1 


100 


Miller-Scalo 


0.354 


0.4 


0.1 


1 




0.354 


1.5 


1 


10 




2.02 


2.3 


10 


100 


Kroupa 1993 


0.579 


0.3 


0.1 


0.5 




0.310 


1.2 


0.5 


1 




0.310 


1.7 


1 


100 


Chabrier 






0.1 


1 




0.238 


1.3 


1 


100 


Kroupa 2001 


0.449 


0.3 


0.1 


0.5 




0.224 


1.3 


0.5 


100 


Baldry & Glazcbrook 


0.323 


0.5 


0.1 


0.5 




0.199 


1.2 


0.5 


100 


Paunchy 


0.315 





0.1 


0.5 




0.194 


0.7 


0.5 


4.0 




0.676 


1.6 


4.0 


100 


Obese 


0.083 


1.35 


0.188 


100 




0.905 


1.65 


1.5 


100 


Extreme top-heavy 


0.136 


0.95 


0.1 


100 



•^ The IMF is written as cW/dlnM = K{M/ Mq)'^, 
for Ml < M < Mu ■ The Chabrier IMF has a lognormal 
form dN/dln{M) oc exp[-(log(M/0.08))2/(2 • 0.69^)] 
below 1 Mq , and is continuous at 1 Mq . In some cases 
the original versions of the IMF continued to lower 
mass, but we have imposed a uniform lower mass limit 
of O.1M0. 



conventional choice is the original Salpeter IMF, with a mass 
range 0.1-100 Mq. However, there is evidence that this con- 
tains too many lo w-mass stars, as sho wn both by local field 
star observations JGould. Bahcall. fc Flvnn 19961 and by es- 
timat es of the dynamical mass in ellipticals (IBell fc de ,Tond 
I2OOII) . The latter authors advocate a "diet" Salpeter IMF, 
which they create by truncating the IMF at the low end 
so that it contains 0.7 times the mass in the current stel- 
lar population, when normalised to the same high-end am- 
plitude. Using a representative estimate of the cosmic star 
formation history from our best fit below, we find that this 
implies 0.788 times as much mass in the IMF before the 
high-mass stars burn away; the indicators of star formation 
rate per mass of stars formed are thus shifted by the in- 
verse of this factor. (This scaling is not very sensitive to 
the history.) These two IMFs and some other forms that we 
consider later in this paper are summarised in Table Q 

To convert between the star formation history and ob- 
served measures such as the stellar mass, UV continuum, 
and - ft'-band light, we use the PEGAS E.2 stellar population 
code JFioc fc Rocca-Volmerangg|l997il .^ This code contains 
a large number of options for the stellar physics. Our con- 
versions assume a close binary mass fraction of 0.05, evo- 
lutionary tracks with stellar winds, the SNII model B of 
Woosley and Weaver (1995), and a constant metallicity of 
Z = 0.02. To assess the sensitivity of our results to this par- 
ticular stellar population code, we have compared them to 



those obtained using the GALEXEV-2003 code of Bruzual 
& Chariot, and find only small differences of ~ 10% at most 
(see also Bruzual & Chariot 2003, where some additional 
comparisons between the codes are presented). 

When constructing the cosmic SFH one needs to con- 
vert from observed luminosity to star formation rate (SFR) . 
Table |5| lists the conversion factors of different star for- 
mation indicators for three different IMFs. The conver- 
sion factors are fboi = Lboi/SFR, fua = Lna/SFR, and 
fuv = Lv/SFR. The latter quantity is nearly constant 
over the wavelength range 1500-2800 A used in various UV 
galaxy surveys, since L^ itself is nearly constant before ex- 
tinction for burst lengths of ^^ 10*-10® yr; hence we choose 
to evaluate fuv at 2000 A without loss of generality. We cal- 
culate these constants with PEGASE, using constant-SFR 
bursts of various duration. For a given SFR, the observed 
luminosity depends on how long that SFR rate persisted be- 
fore the observation. Since surveys are flux limited they are 
biased towards observing galaxies undergoing rapid bursts, 
boosting their luminosities. Therefore, one usually assumes 
a short duration burst when making the conversion to SFR. 
In Table |5| we show the conversion factor for three different 
assumed burst lengths of 10^, 10*, and 10^ years. As shown 
in the table, this uncertainty in burst duration introduces 
uncertainties in the derived SFRs of up to a factor of two. 
In our plots we assume burst lengths of 10* years. We note 
that assuming such globally short burst lengths cannot be 
valid when integrating over the entire galaxy population, as 
is done to get the total SFR at any epoch, since this includes 
galaxies that were too faint to be included in the survey and 
hence had no bias to be observed while bursting. However, 
correcting for this bias with a fully self-consistent approach 
is beyond the scope of this paper. 

Furthermore, the UV flux depends somewhat on metal- 
licity. Using PEGASE we find that for metallicities of Z = 
0.002 and Z = 0.008 the UV output per unit of star forma- 
tion is roughly 25% and 10% larger, respectively, than for 
the solar value of Z = 0.02, for Gigayear bursts. However, 
IPanter et al.l ll2003r found that the mean metallicity of star- 
forming gas has been close to solar for ~6 Gyr into the past, 
suggesting there is not much of an offset in the average UV 
calibration for much of the Universe's history. We estimate 
that metallicity effects can introduce systematic calibration 
uncertainties of ~ 10%, but we will not attempt a more 
detailed treatment until the cosmic metallicity history be- 
comes better known. Finally, a change in the upper cutoff 
mass, from 100 to 150 Mq, for example, would change the 
bolometric and UV output by only a few percent. The SFR 
tracers are produced mostly by stars with M > 5-30 Mq, 
depending on which indicator is used. 

Another required quantity for our study is the mass 
remaining in a stellar population, after the higher-mass stars 
burn away. We use PEGASE to calculate a table of total 
stellar mass versus time after a burst of star formation. We 
include in this total white dwarfs and neutron stars, which 
at late times can be a marginally significant fraction of the 
mass. We then convolve this table by an assumed cosmic star 
formation history /9*/orm(z) to get the remaining present-day 
stellar mass density, i.e. 
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Table 2. Conversion Factors for Star Formation Rate Indicators 



IMF 



lO'^ 10^ lO^yr lO'^ 10* lO^ yr (lO^iergs"! Af ~ Vr) 



Salpeter 


4.31 


6.51 


8.60 


5.00 


7.68 


8.76 


Diet Salpeter 


5.52 


8.33 


11.0 


6.35 


9.75 


11.1 


Kennicutt 


4.94 


8.05 


11.50 


5.78 


9.64 


11.4 



1.22 
1.55 
1.30 



The "mass-lockup" fraction /mass = /5*rem/p*/orm is very 
sensitive to the low-mass form of the IMF and only slightly 
sensitive to the star formation history, for reasonable choices 
of both. For our standard star formation history derived 
below, we find lockup fractions of fmass = 0.71 and 0.63 
for Salpeter and diet Salpeter IMFs respectively. This mass 
is dominated by main-sequence stars with M ^ 0.5 Mq in 
the latter case. However, as noted in the introduction, the 
low-mass form of the IMF sets the calibration between any 
of our luminosity measures and the total stellar mass, but 
it has almost no effect on the relative values of the three 
measures, and therefore on our eventual conclusions. 

The K-h&nA light produced by the current stellar pop- 
ulation is given by a similar equation. 



PK = 



P*form{z) Ik [to - t{z)] dt . 



(2) 



Here Ik is the if-band luminosity per unit stellar mass of 
a stellar population as a function of its age. For a specified 
IMF, ptrem and pk are very tightly correlated since they are 
both dominated by low-mass stars. The if-band luminosity 
is dominated by red giant branch (RGB) and asymptotic 
giant branch (AGB) stars with M « 1.0 Mq. At a lower 
metallicity oi Z = 0.002, we find pK is lower by ~ 5% for 
plausible star formation histories, because of the bluer colour 
of low-metallicity giants. 

We obtain results in J and z bands by a similar 
method. pK as given by Equation|21is the unextincted light, 
but even in if-band there is a small amount of extinc- 
tion. For an extinction estimate, we adopt the "fiducial" 
du st model from the gl o bal lu minosity density modelling 
of iBaldrv fc Glazebrookl ll2003r . This model implies effec- 
tive optical depths in the K, J, and z bands of tk = 0.08, 
Tj = 0.15, and r^ = 0.22. Some idea of the uncertainty 
in these values can be obtained from the three models in 
Figure 6 of Baldry & Glazebrook, along with their fiducial 
model; the dispersion in the values from these four models is 
0.03, 0.05, and 0.07 for the K, J, and z bands, respectively. 

Given the star formation history of the universe pi,{z), 
we can easily compute the resulting bolometric background, 
which we will refer to interchangeably as the extragalactic 
background light (EBL). The background is given by an in- 
tegral over cosmic time: 

CI. c - 

JeBL = -. — / P*form{z) SeBl{z) dt = - — p*formSEBL ■ (3) 

Here the background contribution function Sebl{z) gives 
the total background light energy per unit mass formed into 
stars at a redshift z: 

hoi [tern — t[z)] 



Sebl(z 



■dte 



(4) 



We use PEGASE to calculate the bolometric luminos- 



ity per unit mass hoiit) for a given choice of IMF. For stars 
with lifetimes much shorter than the lookback time to'-t(z), 
the contribution to the integral is simply Eboi{M){l + z)~^, 
where Ei,oi{M) is the total energy released from nucleosyn- 
thesis in a star of initial mass M over the lifetime of the 
star. For less massive stars, the contribution is enhanced 
by the larger scale factor [1 -|- z{tem)]~^ , especially during 
the more luminous, post-main-sequence phase at the end of 
the star's life. At still lower masses, where the stellar life- 
time exceeds the lookback time and the stars do not have 
time to release most of their energy, the contribution then 
diminishes rapidly. Over the entire range of stellar masses, 
Ei,oi{M) is roughly proportional to the mass. Hence, the 
contribution to the background roughly tracks the integrated 
mass in the IMF above the main sequence turnoff. For an 
IMF of Salpeter slope, about half of the total energy that 
the stellar population will give up within the present-day 
age of the universe is emitted within the first 0.13 Gyr, or 
above a turnoff mass of 4.6 Mq. The background light is 
biased to slightly lower masses; for a realistic SFH, half of 
the background comes from stars above ~3Mq. 

The background contribution function Sebl{z) is max- 
imised at a lookback time of about 2.3 Gyr or a redshift 
of 0.2, for a Salpeter or diet Salpeter IMF with a peak 
value of 1.0 x 10^° Lq M" Vr for a diet Salpeter IMF. 
Despite the (1 -I- z)~^ dimming, Sebl from high redshift 
star formation remains above 0.3 times its maximum value 
owing to contr ibutions from long lived, lower mass stars. 
JMadau fc Pozzetti 2000 contains some plots of Sebl(z) for 
Salpeter-like and top-heavy IMFs.) We note that Bruzual & 
Chariot models give very similar results for the bolometric 
luminosity; specifically they are larger by about 4% on aver- 
age, which is well within our errors. The derived bolometric 
background can also be affected by metallicity; it goes up 
by ~ 20% when the metallicity is decreased from solar to 
Z = 0.002, assuming plausible SFHs, though as noted above 
the average metallicity is roughly solar many Gigayears into 
the past. 



3 OBSERVATIONAL MEASURES 

3.1 Total background radiation 

Figure shows measurements of the EBL as a function of 
wavelength. This plot shows absolute measurements of the 
EBL (solid symbols), and estimates resulting from the in- 
tegrated light of galaxies (open symbols). The latter are 
lower limits to the background flux, since faint or low- 
surface-brightness galaxies or exotic (non-stellar) sources of 
energy may also contribute to the background light. The 
EBL is dominated by two main peaks in the optical and in- 
frared, corresponding to direct and dust emission. For now 
let us assume that all of this energy ultimately derives from 
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the radiation of stars. The units we use for both the dif- 
ferential (vJv) and total (Jebl) background radiation are 
"bgu"= nWm^^sr^^ = 10^^ ergs^^ cm~'^ sr~^. A compre- 
hensiv e review of the background is given in lHauser fc Dwekl 
J200ll). 



Integrated gala xy counts in the o ptical /NIR region 
were compiled by iMadau & Pozzettil i200Cl) . who sug- 
gested that the flux had essential ly conv erged. In contrast, 
iBernstein. Freedman. fc Madorj i2002lJl argued that the 
isophotal magnitudes used in these estimates failed to cap- 
ture much of the flux from the faint portions of galaxies, 
and found substantially higher estimates from an "ensemble 
photometry" techni que, as shown by the open circles in the 
figure. iFazio et al.l 112004) provide integrated galaxy counts 
in the NIR from Spitzer, which are not as close to conver- 
gence as those in the optical. 

Direct measurements of the absolute background by 
Bernstein, Freedman, & Madore (2002a, b) give substantially 
higher values. These authors attribute the differences in the 
optical region to two effects: fiux missed by normal galaxy 
count surveys in the outer, diffuse parts of galaxies (par- 
tially recovered by the ensemble photometry technique) and 
high-redshift galaxies that are missed altogether by these 
surveys owing to their low surface brightness. Direct mea- 
surements in the NIR, usin g the space instrum ents DIRBE 
JCambresv et ahlBo Ol: Gori ian et al.ll2000l: fWMe ht & Reese 
I2OOOI: IWrighti 120011) and IRTS iJMatsumoto et alJ I2005D . 
mainly give even higher estimates, particularly the latter. 
The dispersion among these DIRBE points is mainly from 
different ways of analysing the foregrounds, including differ- 
ent ways of subtracting stars and different models for the 
zodiacal light. The zodiacal light subtraction is likely the 
largest source of systematic error in these direct background 
measurements, especially since the shap e of the NIR exces s 
nearly matches that of the zodiacal light JDwek et al.l2005l) . 
It is perhaps encouraging that all of the optical and NIR di- 
rect background estimates from different groups give values 
significantly higher than the direct counts, suggesting that at 
least part of the excess signal is real. To many authors, the 
apparently peaked NIR excess suggests a Lya-dominated 
contr ibution from Pop I II stars before reionisation at z > 10 
('e.g.. lSantos et all2002^ ■ However, this model has some ma- 
.ior problems stemmi ng from its large energy requirements 
iJMadau fc Silkll2005h and observation al constraints on the 
number of NIR-detected Lya emitters llSalvaterra fc Ferraral 
I2OO6I) . 

Turning to the mid- and far-IR region, we show inte- 
grated galaxy counts derived from Spitzer observations at 24 
^m (Papovich ct al. 2004). We also show the inte grated light 
derive d from an analysis at 70 and 160 fj,m by iDole et alJ 
l|2006h . who stacked the observed 24 /im sources and inte- 
grated the resulting light in these other bands. In both cases, 
we show only the fiux resulting from the observed counts, 
rather than their extrapolation of the counts to zero flux, 
which would raise the plotted points by 20-40%. The faint 
counts in these analyses are heavily dependent on estimating 
the detection completeness against the background of unre- 
solved sources, and future work on this issue may well result 
in slightly different faint-end slopes, which would strongly 
affect the extrapolation. 

The absolute background detections in the far-infrared 
(FIR) region come from the DIRBE and FIRAS experiments 



on board COBE . For t he DIRBE data originally obtained 
by iHauser et alj il998l) , different groups have treated the 
photometric calibration and foreground subtraction in dif- 
fcrcnt ways, leading to multiple re sults from the same data 
(Finkbeiner et alJl200d : iLagache et al. 2000: Wright JOoJ. 
At 240 /xm the different results are nearly in agreement, 
but as the wavelength decreases the scatter among them in- 
creases, probably because of the greater complications asso- 
ciated with the higher zodiacal light intensity. These shorter- 
wavelength points can have a strong effect on the estimated 
height and width of the FIR peak. We have represented 
the longer- wavelength FIRAS data by the analytic flt in 
iFixsen et al . ( 1998) , rather than the noisier raw data. 

TeV 7-ray measurements towards blazars can in prin- 
ciple measure the background via the opacity to 2-photon 
scattering. In our judgement, these measurements are not 
yet very reliable; the estimated optical depth is extremely 
sensitive to the assumed shapes of both the intrinsic TeV 
source and the background spectrum itself. Hence at this 
point, we prefer to use direct absolute background detec- 
tions, even though they are uncertain. However, in the MIR 
range only loose upper limits exist as the ga laxy signal is 
too l ow to compete with the zodiacal light IIHauser et alJ 
Il998h . We thus show res ults of optical depth mod els towards 
two blazars presented in A haronian et alj ll2002f) . using the 
formulae of Aharonian (2001J) to convert optical depth to 
MIR background intensity. These arrows are tentative de- 
tections but should probably be treated as upper limits. We 
note that TeV 7-ray absorption estimates should be best 
at constraining sharp, intense features in the background; 
it would be at least somewhat surprising if the NIR peak 
show n by the IRTS data is consistent with existing TeV 
data i Aharonian et al]|2005^ . To strengthen the case for a 
lower EBL in the MIR region, we show the results of fluctua- 
tion an alyses by Kashlinskv et al. ( 1996a j, Kashl inskv et alJ 
l|l9963), and Kashlinskv & Odenwald (2000), which at- 
tempt to place limits on the galaxy contribution from the 
power spectrum of the background rather than its absolute 
intensity. 

As mentioned above, and as is apparent in the galaxy 
counts in Figure the galaxy contribution to the extra- 
galactic background is expected to ha ve two broad humps . 
A model of the background radiation JPrimack et alj|2005h 
is shown by the dotted line; while this clearly underesti- 
mates the background intensity especially in the FIR, it in- 
dicates the expected shapes. For our work, we require bet- 
ter estimates of the total intensity. We show three traces 
through this plot, representing what we consider the min- 
imum (dashed line), maximum (dot-dashed line), and best 
guess (solid line). These curves are simply constructed to 
be smooth curves without a particular functional form, but 
their shapes do reflect global galaxy emission models. The 
"max" trace uses the background measurements in the opti- 
cal and NIR and the MIR gamma-ray limits, while the "min" 
trace follows the lower values implied by the galaxy counts. 
The "mid" trace adopts a compromise between the higher 
galaxy counts and the lower range of absolute background 
estimates. In the FIR, the differences between these curves 
are unrelated to those in the optical, and mainly reflect the 
assumed width of the FIR peak. The total background fluxes 
for these traces are 50, 77, and 129 bgu. If we consider only 
the optical/NIR portion below 10 fj,m, we get 22, 36, and 70 
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Figure 1. The observed extragalactic background light. Sources for observational points and limits are given in the legend. Filled 
symbols represent absolute measurements of the background light, while empty symbols represent integrated galaxy counts. The filled 
circles for the Bernst ein et al. (2002a, b) absolute background incorpor ate the changes in their erratum. We show the analytic fit to the 
FIR measurement of lFixsen et al.l il998r . For the gamma-ray limits of lAharonian et al.l i2002^ we show squares corresponding to their 
tentative detections of t = 1.5 at 5.5 TeV for Mkn 501 and t = 5.5 at 4.0 TeV for H 1426+428. We con vert these limits to effective 
IR wavelength and background inte nsities using the approximations of lAharonianI i200U) . Upper limits from lKashlinskv et al.l il99(ialhh 
and IKashlins kv fc Odenw ald f2000|) use the power spectrum to estimate the total extragalactic contribution. The dashed, solid, and 
dot-dashed black curv es sh ow th e traces th rough this data that are discussed in the text. The dotted curve shows results from the 
semi-analytic model of lPrimack et al.l J2005|). 



bgu, while if we consider only the FIR portion above 10 /ini, 
we get 28, 41, and 59 bgu. 

Our EBL estimates range to substantially higher val- 
ues than the best estima te of 55 ± 20 bgu determined by 
iMadau fc Pozzettil i2000h. which was lat er revised and in- 
creased to 60 bgu bv lMadau et a l. (2001). Later estimates, 
which include 60-93 bgu by jGisport, Lagache, & Puget 
J2000h. 45-170 bgu with a preferred value of 100 

bv iHauser fc Dwelj J200jl. a nd 100 ± 20 bgu by 

iBernstein. Freedman. fc Madorg ll2002lJ) . are more in line 
with our estimate. The crucial issues contributing to this 
uncertainty are not statistical error but systematics like the 
caUbration of the FIR measurements, the zodiacal Ught sub- 
traction in the optical and NIR background measurements, 
the treatment of the falloff from the optical and FIR peaks 



to the MIR region, and simp ly which da t a poin ts to include. 
Illustrating the latter issue, iDole et al] (|2003) recently dis- 
cussed the EBL in the context of their FIR galaxy stacking 
analysis. They rely on limits from TeV 7-ray absorption and 
IR fluctuation analyses to set the upper limits on the flux, 
whereas we have preferred absolute background detections. 
The upper and lower boundaries of their allowed region are 
shown in their Figure 12. Using these traces and adopting 
a higher-order integration scheme than used in their deriva- 
tion, we flnd that the Jebl is in the range of 52-77 bgu, 
spanning approximately the lower half of the allowed range 
we flnd here. 

Below, we will need to estimate the probability (or like- 
lihood) distribution of Jebl- It is clear that this distribution 
is far from Gaussian, as the galaxy counts provide a hard 
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lower limit to the background in certain portions of the spec- 
trum, but the correlated and conflicting data points mean 
that it is hard to estimate the likelihood function in an ob- 
jective way. Our approach is to split the wavelength regime 
in two at 10 fim, since the values in the optical/NlR and 
FIR are essentially unrelated. In each of these ranges, we 
take the likelihood to have a constant (top-hat) distribution 
in log space between our minimum and maximum traces. 
We then convolve the two distributions together to describe 
the sum of the flux in the two wavelength regions. 

Of course, active galactic nuclei (AGN) also contribute 
to the background energy. This energy is emitted mainly in 
a UV/optical peak, a roughly power-law X-ray tail, and an 
MIR/FIR peak from dust-processed radiation. However, the 
AGN contribution to Jebl appears to be mu ch smaller than 
that o f star-forming galaxies. For example, iHopkins et alJ 
i2006 h present a detailed model of the AGN and supermas- 
sive black hole populations that fits (partly by construction) 
quasar luminosity function data over a wide range of wave- 
lengths and redshifts. Their Figure 22b plots the growth of 
the average black hole mass density with redshift, and we 
can compute the corresponding contribution to the bolo- 
metric background radiation using their assumed radiative 
efficiency of tr = 0.1 (essentially using the argument of 
ISoltan|[l983) . We obtain Jagn = l-9loJ bgu, where the 
uncertainty is scaled from the uncertainty in their predic- 
tion of the z = black hole mass density, pbh ~ 2.9ljg x 
10^ Mq Mpc~^. Note that th is argument does not assume 
that the lHopkins et alJ (120061) scenario is physically correct, 
just that it fits the observed luminosity functions within 
their uncertainties and, therefore, reproduces the observed 
total emission. Furthermore, the mode l's predicted pbh 
agrees with observational estimates fe.g.. lAller fc Richstonel 
l2002l : iMarconi et al.ll2004) , and there is no room for a large 
amount of additional AGN contribution to Jebl without 
overproducing these estimates. A significant fraction of this 
bolometric AGN background appears in X-ray wavelengths 
and therefore does not add to the UV/Optical/IR back- 
ground considered here. A s another check of Jagn, the 
model of ISilva et alJ i2004) . which is based on AGN counts 
in X-ray and IR bands an d has recently received further ob- 
servational support from iTreister et ahl (|^006), implies the 
AGN contribution in all wavebands is ~ 1.4 bgu, and that in 
the UV/Optical/IR background it is 1.1 bgu (Silva, private 
communication). We conclude that the AGN contribution 
cannot be much larger than 4%, and is probably somewhat 
lower, and we ignore it henceforth. 

The issue of whether or not the integrated galaxy counts 
match the total EBL has implications for this paper beyond 
just a change in the total bolometric background. If flux 
from sources besides normal galaxies are truly required to 
explain the extragalactic light, and the source cannot be 
identified, then there is not necessarily a relationship be- 
tween the total EBL and the energetic output from stars! 
However, the need for additional sources of energy is not at 
all clear at present. The systematic uncertainties, particu- 
larly in the galaxy photometry and zodiacal light subtrac- 
tion, suggests that a curve splitting the difference between 
the raw counts and measured total background (as in our 
"mid" model) represents the most conservative assumption. 
Also, if the excess over total normal-galaxy counts should 
turn out to be real, a candidate for explaining the excess 



in the NIR, where it appears most significant, exists in the 
form of Population III stars. According to models of these 
sources, they would contribute almost all of their energy in 
the wavelength range 1-5 /^m, leaving the background at 
other wavelengths essentially unchanged. This would imply 
that our lower limits on the optical/NIR background (which 
ignore the NIR excess) and our traces of the FIR background 
are still valid. 



3.2 Star Formation History 

The cosmic star formation history is related to both the 
EBL and the stellar mass observable today. Observational 
estimates of the SFH at various redshifts are shown in Fig- 
ure|21 This plot uses a combination of IR (stars), UV contin- 
uum (filled symbols), and emission-line surveys (open sym- 
bols) . We have tried to use the deepest and broadest surveys 
within each type, so some of the pioneering SFH surveys 
have been omitted. The plot also shows the sub-mm points 
from iBarger. Cowie. fc Richarda ll200d . crosses) , including 
both their actual source detections and their extrapolation 
to account for the entire sub-mm background, as well as 
the estim ate of the total SFR in the damped Lya popula- 
tion from I Wolfe et al.l (l2003ll . We do not use the sub-mm 
or damped Lya points in the fit. We omit estimates from 
radio and X-ray surveys, because the energetic calibrators 
are even more uncertain than those in the UV and FIR. 

For the UV continuum m easurements, we have com- 
bined resu lt s at 2 800 A from IConnollv et al.l lll997ll and 
Wolf et al.' ("2003^, and re sults a t 1500-170 A fi-pm 
Steide l et al. ( 1999), Gabasch^elai] ll2004'l. Giavalisco et al. 
B004), ISchiminovich ct al. (2005|), ISawicki fc Thompson 
[20051) , and iBouwens et aL (■200&1 . For the IR measurements 



we inc lude the local IRAS survey of|Yui^^ed^^^^ondpn 
20W ) and the Spitzer 24 /im survey of lPerez-Gonzalez et al. 
20051). For emission-line measurements, we restrict our- 



selves to Ha because metal lines like [O II] might be biased 
b y evolution of th e metallicity and we use the values listed 
in lHopkinj i2004l) from t h e surveys o f^allcgo et al. ( 1995* 



Trossc & Maddox 



X 



(J200^ from t h e surveys o f|Gc 
addox (1998'). Y an et alJ IU9 



19981) .ISullivanet al 



20001. Hopkins et al (2000h . iTresse et al.1 J2002D . and 
iPerez-G onzalez et alJ i 2003l) . We use only the Ha estimate 
from lSuUivan et al.l i200Ci) because the UV continuum mag- 
nitudes used there have been brought into question by more 
recent results from GALEX dSchiminovich et al]l2005^ . 

We convert the calibration factors used by the vari- 
ous authors to the uniform values given in Table |5| using 
a diet Salpeter IMF with 10* yr bursts. We have converted 
all the star formation rate densities to our standard cosmol- 
ogy, since they scale with the Hubble constant, H{z) = d/a. 
This cosmology dependence factors out when integrating 
over time to get pi,form , although p*rem can shift by a percent 
or so owing to the dependence of fmass on cosmic timescales. 

The UV surveys, in general, do not probe deep enough 
in luminosity to make the integrated star formation rate con- 
verge. Most authors calculate the total by fitting a Schechter 
form to the luminosity function and extrapolating either to 
zero flux or to an arbitrarily chosen lower limit. For consis- 
tency, we adjust these reported values to a uniform extrap- 
olation down to O.IL* where L, is the Schechter turnover 
luminosity reported for each survey. Extrapolating down to 
O.OIL, instead would cause the values to increase by 20- 
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Figure 2. Star formation rate per unit conioving volume as a function of redshift as derived from observations. The observational points, 
which are labelled by source, have been rescaled to our chosen "737" cosmology, and to a common set of luminosity and completeness 
calibrations using the "diet" Salpeter IMF as discussed in the text. We use uniform dust correction factors for the optical and UV points 
(see the text for details). The thick solid line is our best fit to the data. The thin solid and thin dashed curves indicate the la envelope 
of the fits in our standard and "random-calibration-error" samples of fits, respectively. 



45%. In many cases ours is the approach used by the orig- 
inal authors; in the others the change in the extrapolation 
has a rather minor effect compared to the scatter between 
surveys. In contrast to the UV, the Ha surveys are deep 
enough that they require no extrapolation. 

To make the UV continuum extinction corrections rea- 
sonably uniform, we use the correction formula of Calzctti 
il999l). assuming the m ean value of E{B — V) = 0.15 given 
bv lSteidel et alJ (|l993)- This yields correction factors of 3.1 
and 4.7 at 2800 and 1500 A, respe ctively. The one excep- 
tion is the highest-redshift point of iBouwens et all i200£tl . 
for which we have used the correction factor of 2.0 estimated 
by the authors from the mean UV slope in the sample. At 
lower redshifts, some authors have used this UV slope cor- 
rection method to obtain higher extinc tion values than we 
have assumed (e.g., a factor 7 in both ISchiminovich et alJ 
l2005l and Giavalisco ct al. 2004i ). 

For the Ha measurements, we used the uniform 



luminosity-dependent extinction calibration of iHopkina 
(2004). The two FIR surveys are not "transmission- 
corrected" for the amount of UV light that leaks out directly 
from the galaxies, but this is likely a significant fraction. In 
fact, at low redshi fts the UV and FIR luminosity estimates 
are roughly equal IIMartin et al. 2005') . We have crudely as- 
sumed a correction factor of 2.0 for 2 < 1 and 1.2 for z > 1. 
All of these extinction and transmission estimates are clearly 
highly uncertain, and among the largest of the calibration 
errors in the observed SFH. 

In Figure 13 we replot the SFH such that one can eas- 
ily integrate under the data points by eye to get the total 
mass of stars formed: dpi, /da versus scale factor a. This plot 
makes it straightforward to see the relative importance of 
star formation at different epochs in determining the stellar 
content today. Clearly, most of the uncertainty in the total 
stars formed comes from the redshift range 1 < z < 3. While 
there is still a surprisingly high dispersion in low-redshift 
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Figure 3. Star formation rate per unit scale factor versus scale factor. This is similar to Figure|^but with rescaled axes. The vertical 
axis here is the derivative of the mass of stars formed with respect to the scale factor, dp^^f^j-^/da = (1 + z)pf^fgj.^/H{z). In this form the 
data points are independent of cosmology, and the total star formation is the area under the data points. One can easily see the relative 
importance of different redshift ranges. 



{z < 0.5) estimates of the star formation rate, the figure 
shows it is of httle significance to the total stellar mass. 

In these plots, it is unclear how well individual ob- 
servations capture the whole census of star formation 
at any redshift. At high redshift, several populations of 
sources have been suggested as major contributors to 
the total star for mation, including diffuse red galaxies 
iWebb et alj l200ai . damp ed Lyo: absorber s i Wolfe et alJ 
l2003l) . L yg emitters (e.g.. iTapke n al.ll20 0(Jl. and sub-mm 
galaxies (IBarger. Cowie. fc Richards! 1200(1 : IChapman et alJ 
12005). For all of these categories, the extent of overlap with 
the Lyman-break galaxies, for example, is still unclear. The 
mere presence of one population in a second population 
may not suffice; for example, the large star formation rates 
found for the bright sub-mm sources are e ssentially invisi- 
ble to the SFR estimates made in the UV iChapman et alJ 
|2003)- Hence, it is possible that the SFR estimates obtained 
from these different populations should be added together, 
rather than plotted together as different measures of the 



same quantity. This could raise the true SFR by perhaps a 
factor of 2 above the individual points. 

To compare this large spread of data points to other 
measures of star formation, it is useful to obtain analytic 
fits to the star formation history. Here we introduce a new 
parametric form: we set 



P*form{z) = po/(l +Pia ^^Y^ 

where a = 1/(1 -I- z) is the scale factor. Then 






(1 +pia-P2)P3H 



(5) 



(6) 



This functional form has the advantages of easily relating 
p-iform and pi: form, aud of making the dependence on cosmol- 
ogy explicit. This parameterisation is a particularly natural 
choice for Figure |21 since dpt^jorm/da = p^jorml \<^ii(^z)\ is a 
simple analytic function. 

We use a x'^ technique to fit the data. The reported er- 
rors for the different points vary widely, depending mostly on 
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sample size; for those with the smallest errors, the system- 
atic errors in calibration and extrapolation probably swamp 
the formal errors, which are derived from Poisson statistics 
alone. Rather than trying to revise the error derivation of 
every literature source, we simply add an error of 0.1 dex 
in quadrature to the errors on the points before fitting the 
data, which restores the best-fit x^ to an acceptable value 
of 44 for 42 degrees of freedom. The best-fit parameters 
are po = 9.0 x 10* M© Mpc-^ pi = 0.075, p2 = 3.7, and 
P3 = 0.84. 

To test for a possible dependence of our results on th e 
parametric form, we also use the form of ICole et al.l i2001l) : 



p-kjor 



(a -I- hz)h 
1 + [zlcY 



(7) 



Here the dependence on the present-day Hubble constant 
is factored out, although not that of the other cosmolog- 
ical parameters. This form has the same number of free 
parameters as the previous one and has been widely used 
in the literature. Fitting the data as before, our best fit 
is a = 0.0103 Mgyr-^Mpc"^, h = 0.088 Mq yr"^ Mpc"^, 
c = 2.4, and d = 2.8 with x^ = 44.7 for 42 degrees of free- 
dom. The curves produced by these two parameterisations 
are nearly indistinguishable from each other. 

We want to probe not just the mean value but the full 
allowed range of star formation histories. To do this, we con- 
struct an ensemble of parametric fits, using a Monte Carlo 
Markov chain and assuming uniform priors on the logarithm 
of the parameters. The ensemble evenly samples the prob- 
ability distribution of the fits. We can then estimate the 
allowed ranges of p*/orm(z) at any z from the ensemble of 
fits. (A simple grid search in parameter space was also tried, 
and produced essentially identical results.) The typical for- 
mal 1(7 error in p*/orm is only 10-15% for 2: < 5, as shown 
in Figure 121 The width s of these envelopes are quite similar 
to those of Hopkins fc Beacom ( 2006) , despite the difi'erent 
data sample, sampling procedure, and parametric form. 

However, it is likely that this procedure underestimates 
the true error s in the determination of P -khrm,- For example, 
the best fit in lHopkins fc Beacoml i2006r appears to lie well 
above our error envelope for 2 < z < 4, owing to a different 
selection of data points, different assumptions about extinc- 
tion, and a larger extrapolation of the luminosity function. 
While the large number of data points with small Poisson 
errors result in a small formal error in our fit, this masks 
the larger uncertainties in calibrating /9*/orm, which applies 
equally to all the points. Even with a fixed IMF, there are 
substantial uncertainties in the mean opacity at a given red- 
shift and in the calibration of luminosities owing to the un- 
certain star formation histories of individual galaxies (see 
§2). The systematic errors at different redshifts are undoubt- 
edly correlated, although there is speculation that the mean 
opacities at « ~ 0, « ~ 2, and z > 4 are substantially differ- 
ent from one another. 

Therefore, we produce an alternate grid of fits in which 
we randomly and smoothly adjust the p^form data values 
as a function of redshift. We first define three random vari- 
ables A{Q), ^(1.5), and A(5), and set each to have normal 
distributions of dispersion ua ~ 0.5. We then multiply the 
data values and associated errors by e '^', where A(z) is 
the quadratic function in z through the three points ^4(0), 
A{1.5), and ^(5), and repeat the previous Monte Carlo sam- 



pling procedure. Not surprisingly, this procedure results in a 
much broader range of acceptable parameters, and a corre- 
spondingly wider envelope of p^foTm at any given redshift, as 
shown in Figure |5| Given the large uncertainties discussed 
above, oa = 0.5 is probably reasonable. The I-ct dispersion 
in the SFH is then about 50% or 0.2 dex for z < 4. 



3.3 Stellar mass and NIR luminosity density 

The third observable we need to consider is the stellar mass 
density. We are primarily concerned with the value of this 
quantity today but we shall also consider the evolution of 
stellar mass with cosmic time. NIR light is a useful proxy 
for stellar mass, since plausible values of M/L vary only by 
a factor of a few and this variation can be acc ounted for 
using simple colour-based estimates JBell fc de J one 200 ij). 
Like the star formation rate, the stellar mass and luminosity 
densities determined from observations are proportional to 
H{z), and we have translated the measurements below to 
our assumed cosmology. 

We shall use five studies of NIR light and stellar mass 
at low redshifts. All the studies in the J and K bands use 
2MASS magnitudes, combined with various redshift surveys. 
The results of these lo cal surveys are given in Table |21 

ICole et aO i200lf) use the 2dF redshift survey in the 
southern sky, combined with the second incremental release 
of 2MASS, to obtain J, K, and stellar mass densities. They 
use Kron aperture magnitudes, calibrating the preliminary 
2MASS catalogue against deeper _ft'-band data and correct- 
ing their estimate of p*rem for missing light outside the Kron 
radius. However, they do not similarly correct their reported 
estimate of pk, so we have scaled up their Kron magnitude 
based estimate by 0.11 magnitudes to get the total light. To 
determine the stellar mass density they estimate the mass- 
to-light ratio galaxy by galaxy using the observed optical- 
to-near-infrared colours to constrain the effects of the star 
formation history and metallicity in each galaxy. They check 
for the effect of missing low-surface-brightness galaxies in 
the relatively shallow 2MASS sample and conclude that that 
any such effect should be small. 

iBeU et alJ ll2003t) use the Sloan Digital Sky Survey 
(SDSS) Early Data Release in the northern sky, combined 
with the final 2MASS survey, to obtain A'-band and stellar 
mass densities. They also use Kron magnitudes but scale 
up the luminosities of the early-type galaxies by 0.1 mag to 
account for light outside the Kron aperture. They estimate 
their M/L ratios by means of optical galaxy colours; they ne- 
glect the effect of extinction entirely since it affects the mass 
estimate both through the flux and the colour, and these 
should nearly cancel out JBell fc de Jonell200lh . They go on 
to suggest that the 2MASS catalogue may miss some low- 
surface-brightness galaxies and estimate the magnitude of 
this effect in three ways. First, estimating the ii'-band light 
density using the r-band luminosity density, which should be 
complete at these magnitudes, and the global r — K colour 
for their sample, they find pK increases by 12%. Second, us- 
ing their <;- band-selected galaxy sample and synthesising the 
K magnitudes of galaxies missing from the ii'-band sample 
from their optical magnitudes, they find pK increases by as 
much as 30%. This g-band estimate might be biased since 
when they obtain the stellar mass density using the com- 
plete g-band optical sample it is only 4% higher than their 
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estimate using their uncorrected K-hand magnitude limited 
sample. Hence, in this paper we use their direct estimates of 
mas s and light from the is'-band sample. 

IPanter et al.l f2004) analyse galaxy spectra in the first 
data release of the SDSS at 20 A resolution. They obtain 
total fluxes by normalising the spectroscopic fibre magni- 
tudes to photometric _R-band Petrosian magnitudes. (We 
note that Petrosian magnitudes may underestimate the flux 
from early-type galaxies by about 20% , although they are 
probably very accurate f or disk galaxies iStrauss et alj2002l : 
iGraham fc Driven |2005^ . They then fit galaxy models to 
compressed versions of the spectra to obtain the local stellar 
mas s density. 

lEke et alJ (|2005|) use an updated version of the 2dF red- 
shift survey together with the updated 2MASS c atalogue to 
produce a larger sample than lCole et al.l (120011) . obtaining 
J-band, isT-band, and stellar mass densities. Eke et al. cal- 
culate galactic stellar masses using star formation histories 
that include bursts. When tested against mock galaxy cat- 
alogues, this reduces the average stellar masses by 20-30%. 
Whether or not this results in a smaller bias, however, de- 
pends on the validity of the star formation histories in the 
mock catalogues. Eke et al. show that the mass function ob- 
tained from their optically-selected sample differs system- 
atically from that obtained using the J-band selected sam- 
ple, which they cite as possible support f or the bias agains t 
low-surface-brightness galaxies found bv iBell et alj i2003t) . 
However, the difference starts starts to occur where the com- 
puted biases in the optically-based and NIR-based mass es- 
timates from the mock catalogues start to diverge, making 
its significa nce unclear. 

Finallv. lJones et al.l J200d) combine 2MASS with a pre- 
liminary catalogue from the wide-angle 6dF redshift sur- 
vey and determine J-band and K-hand luminosity densi- 
ties. They use "total" (extrapolated) 2MASS magnitudes to 
obtain the galaxy luminosities. Since they do not explicitly 
state a error estimate, we read it off their Figure 15. They 
also argue that their sample is robust to surface brightness 
effects on the total luminosity density. 

For our purposes we must compare the different esti- 
mates using the same assumptions. First, we compute the 
luminosity densities for each estimate using the same value 
for the sun's absolute magnitude in each band, obtaining the 
densities in Table|3 Since the luminosity densities were com- 
puted in a similar way by each group, we first average the 
luminosity density estimates after weighting by the inverse 
of its variance, and take the harmonic sum of the variances 
to get the statistical error. We then add an overall system- 
atic error of 15% (cf. iBell et al.l l2003). which dominates the 
error in all cases. 

To compute the stellar mass densities, we first convert 
all the mass estimates to a diet Salpeter IMF, assuming 0.7 
times as much mass as for a Salpeter IMF and 1.74 times as 
much mass as for a Kennicutt IMF. We then take a straight 
average of the four estimates, since for each data point the 
systematic error outweighs the statistical error and we do 
not know which approach is best. We then a ssign an overall 
systematic error of 20% fcf. iBell et alJl2003r . which may be 
reasonable as the varied approaches here yield a dispersion of 
only 10%. Our final estimate of (3.50±0.70) x 10* Mq Mpc"^ 
for a diet Salpeter IMF implies f2*/i = (1.8 ± 0.4) x 10"^ 
and F* = n^/^t = 0.063 ± 0.016. This estimate is lower 



by more than a factor of two than the earlier best estimate 
of iFukugita. Hoean. fc Peeblea il998ll using _B-band data, 
once we account for the different IMF. 

Recently, attempts have been made to estimate the stel- 
lar density at higher redshift as well, based on faint galaxy 
catalogues with multi-band photometry. This is even more 
difficult than estimating the high-redshift star formation 
rate density, both because surface brightness dimming inter- 
feres more with regions of less intense emission and because 
their estimates must be based on the rest-frame optical light 
rather than rest K-hand, which increases the uncertainty in 
the mass. In Figure |1| we show a number of observational 
estimates for p*(z:). Most of these points were originally de- 
rived for a Salpeter IMF. We rescale them to a diet Salpeter 
IMF using the ratio of population masses calculated from 
our best-fit star formation history; we find this to be a factor 
of 0.70-0.74. The different studies use a variety of assump- 
tions about stellar populations, extinction, and methods for 
calculating the mean M/L ratios. In most cases we plot the 
error bars given by the authors, w hich are usually purely 
from Poisson statistics. However, for lPickinson et al.l (120031) 
we plot error bars based on their dominant uncertainty, mod- 
elling the stella r population. The largest estimates of p* at 
high z are from lGwvn fc Hartwicld i2005r but it is not clear 
whether or not these higher estimates arise from their use 
of deeper UDF data. 

The systematic errors in the measurement of the to- 
tal stellar mass and light deserve some further discussion. 
Extrapolation beyond the faint or bright end of the ob- 
served luminosity function is not a significant facto r in the 
surveys of the local uni verse discussed above (C ole et alJ 
I2OOI1 IPanter et al.M2004^ . These surveys make corrections 
for flux outside the apertures and these apertures are quite 
wide to begin with, e.g. ~ 20-40 kpc in radius for a mas- 
sive disk galaxy, but these corrections do not take into ac- 
count the power-law halos that commonly surrou nd galaxies 
JMouhcine et a l. 2005: Guhathakurta et al.l2005^ . However, 
usi ng the mass fraction and den sity profile suggested for M31 
bv lGuhathakurta et alJ (l2005l) . it appears the extra correc- 
tion would only amount to 2%. 

Free stars not clearly bound to any single galaxy 
could be another so urce of error. Observati ons of red stars 
JDurrell et alJl2002h and planetary nebulae (lOkamura et alJ 
I2OO2I) suggest free stars make up 10-20% of all the stars in 
the Virgo cluster. However, as clusters contain only ab out 
2% of all the stars in the universe dEke et alj l2005l) . it 
by itself is not a significant correction. In simulations of 
galaxy formation presented in a related paper (Fardal et 
al. 2006, in prep), free stars are found preferentially in and 
around galaxy clusters owing to the intense tidal interac- 
tions present there, so once again the global fraction of free 
stars should be well below that estimated for Virgo and be 
negligible for our purposes here. 

Systematic effects from low-surface-brightness galaxies 
unrepresented in the 2MASS catalogue are potentially a 
larger problem; the various estimates of their significance in 
the discussion above range from to 30%. We will assume 
that this is not a significant bias, but the reader should keep 
in mind that the true luminosity and mass densities could 
potentially be raised by as much as ~30% by this effect. 

Finally, the calculation of M/L ratios of individual 
galaxies, something which is required in all these papers. 



© 0000 RAS, MNRAS 000, 000-000 



12 Mark A. Fardal et al. 



Table 3. Local Stellar Mass and Light (for diet Salpeter IMF) 



1.93 ±0.20 Bell et al. (2003^ 



P'j 



weighted avg 



2.11 ±0.31 
2.43 ±0.08 
2.22 ±0.15 
2.37 ±0.36 



Cole et al. ('2001') 

Eke et al. ('2005') 

Jones et al. (2006') 

(assuming 15% systematic error) 



Pk 



weighted avg 



4.24 ± 0.64 
4.08 ± 0.09 
4.66 ±0.15 
4.30 ± 0.37 
4.24 ±0.64 



Cole et al. (2001') 

Bell et al. (2003') 

Eke et al. (2005') 

Jones et al. (2006') 

(assuming 15% systematic error) 



unweighted avg 



3.94 ±0.58 
3.71 ±0.07 
3.26 ±0.11 
3.08 ±0.14 
3.50 ±0.70 



Cole et al. (2001') 

Bell et al. (2003') 

Pan ter et al. (2004') 

Ek e et al. (2005') 

(assuming 25% systematic error) 



"Using M0(^) = 4.53 (Bell et al. 2003') 



"Using M, 



&{J) 



■ 3.33 ("Worthcv 1994') 



=Using M0(^) = 3.70 (Worthev 1994") 



is uncertain owing to many factors including uncertainties 
in: the metallicity, the star formation history, the extinc- 
tion, the stellar tracks, and the population synthesis meth- 
ods themselves. The derived distribution of the metallicity 
may also have a substantial effect, particularly on the K- 
band hght from RGB and AGB stars. The latter stars are 
quite difficult to simulate, and regardless of the metallicity 
they are a major source of error in the stellar population 
calcul ation. To indicate t he dispersi on, the differe nt meth- 
ods oflCole et al.l (l200ll) . iBeTl et alj (i2003) . and Eke et alJ 
llooi) give M/Lk estimates of 0.88, 0.92, and 0.62 respec- 
tively, and our combined sample gives a net result of 0.83. 
This range of estimates is in accord with the '^ 25 % system- 
atic er rors in M/L (for a fixed fMF) estimated bv lBell et alJ 
l|2003h . Observationally testing the distributions of star for- 
mation, extinction, and metallicity derived in local surveys 
will help control the systematic errors involved in calculating 
the stellar mass. 



4 COMPARISON OF STAR FORMATION 
MEASURES 

■We now turn to the central question of this paper: are the di- 
rect estimates of the cosmic star formation history, the mea- 
sured bolometric background intensity, and the local density 
of K-hanA light mutually compatible, given the assumption 
of a universal IMF with a Salpeter-like shape above 0.9 A/q? 
The curves in Figure |1] show the stellar mass density as a 
function of redshift computed from the star formation his- 
tory traces shown by the corresponding curves in Figures |21 
and|21 The lowest curve, which reduces the standard normal- 
isation of the SFR estimates by about a factor of two (see 
§3.2), roughly follows the median of the data points, while 
the other curves are systematically higher than the data. 
■We have included the slight variation with redshift of the 
recycled mass fraction in this calculation. Although many 
factors could contribute to systematic errors, it is of partic- 
ular note that a substantial amount of stellar mass could be 



missing from the high-redshift observational surveys owing 
to surface brightness effects. The tension between the SFR 
and stellar mass density estimates provides circumstantial 
evidence for such effects. If lower surface brightness galaxies 
are indeed missed in the high redshift counts, it would help 
resolve the discrepancy between the absolute background 
light estimates and the integrated count estimates, suggest- 
ing that the higher, absolute estimates are more accurate. 

For the remainder of this section, we will focus on z = 0, 
since the observational estimates of the stellar mass density 
are more secure in the local universe and, even more im- 
portant, we can bring in the bolometric background (which 
is only measured at 2; = 0) as an additional observational 
constraint. 'We will mainly focus on the local K-hanA hght 
density pK rather than the stellar mass density p*rem , since 
it is the directly observable quantity. We integrate each SFH 
fit to obtain the predicted final stellar mass p*rem and K- 
band luminosity density pK as in Equation^ and the EBL, 
Jebl, as in Equation (The integrations cover the range 
2; < 15 in each case.) We include the extinction correction 
e^^^ discussed in S|5| to make pK an observed luminosity 
density. 

According to our direct SFR fits (assuming no calibra- 
tion error), the mean p^rem ~ (5.3±0.3) x 10** Mq Mpc~^, or 
-F* = 0.095 ±0.005. The corresponding predicted extinction- 
corrected PK = (6.9 ± 0.4) X lO**L0j^Mpc"^. (We obtain 
the same results whethe r we use ou r new parametric form 
in Equation ini or that of lCole et al] (|20o3) in Equation |7|) 
Including an extra calibration error with a a = 0.5, we ob- 
tain pK = (6.9 ± 2.2) X 10* Lqk Mpc"^. Our adopted value 
for the observed density is (4.2 ± 0.6) x 10** Lq^^ Mpc~^. 
Thus the integrated SFH should produce more mass and 
light than that directly observed, by a factor of 1.5 if we use 
the quoted mass estimates, or by a factor 1.6 if we use the 
Tf-band luminosity density. The probability of agreement is 
< 2% without random calibration errors, but rises to 10% 
and 32% for random calibration error of oa = 0.2 and 0.5 
respectively. Thus, while the offset could indicate that lo- 
cal surveys miss significant amounts of stars, or that some 
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Figure 4. Stellar mass as a function of z scaled to a diet Salpeter IMF (at 2 = this scaling factor is 0.70 and it slowly increases 
to 0.74 at z = 5). Sources of the data points are as indicated. The heavy solid line shows the result of integrating our best fit through 
Figure |3] including the appropriate recycling factor. The light solid and dashed lines show the Icr boundaries of our "normal" and 
"randomly-recalibrated" ensembles of fits, respectively. 



other assumption (such as our IMF) may be wrong, calibra- 
tion uncertainty in the SFH seems a sufficient explanation. 
For example, the curve could be made to match the plots 
by changing the estimates of the UV extinction. Our best 
SFR fits imply M/Lk = 0.7 1e^^ = 0.7 7 . whi c h is within 
the ra nge o f values found b y ICole et al] i200lf) , iBell et all 
J2003D . and lEke et al] J200. 



The appearance of a discrepancy is found by other au- 
thors as well, although a number of authors apparently sim- 
ply scale the integrated SFH to match the z = Q stellar 
mass, without mentioning the scaling factor! Our integra- 
ti ons of p-^rem are in reasona ble agreement with the results 
of iHopkins &: BeacomI ll200(Jl . who found about a 10% un- 
certainty in the integrated stellar mass from their star for- 
mation fits and a value exceeding that expected from direct 
observations by nearly a factor of 2 (for comparisons us- 
ing the same IMF). The discrepancy implied by their final 
mass is even larger than in our case, because their much 



larger SFR in the range 2 < 2 < 4 contributes a significant 
amount of mass. 

However, we can also include the EBL in the compar- 
isons. Our adopted allowed range of Jebl is 50 to 129 bgu, 
with a most likely value of 77 bgu. Our derived value from 
the SFH fits is (63 ± 4) bgu without including random cali- 
bration error, and (62 ± 20) bgu with aA = 0.5. 

After allowing for fairly large systematic uncertainties, 
it appears that the estimated cosmic star formation his- 
tory can be made consistent, individually, with measure- 
ments of the local K-ha.nA light density and the bolometric 
background intensity. One might, therefore, be tempted to 
conclude that these three probes of the cosmic star forma- 
tion are mutually consistent. Figure [5] shows that this is not 
the case. In the upper left panel, we show acceptable fits 
in the space of pK vs. Jebl, assuming no systematic errors 
in the calibration of the SFRs. The vertical lines represent 
the mean and 2cr errors in the local A'-band luminosity den- 
sity from the combined observations, as discussed in § 13.31 
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Figure 5. The combined range of X-band liglit and EBL allowed 
by the range of SFHs. In each panel, the horizontal lines denote 
the minimum, best guess, and maximum background intensities. 
The vertical lines denote the best fit and la error bars on the 
ii'-band luminosity density. The vertically extended rounded con- 
tours show the joint 67%, 95%, and 99.8% joint confidence inter- 
vals on these two observed parameters. The diagonally extended 
contours show the relation of the two quantities in the ensem- 
ble of SFH fits, for the following assumptions; (a) No freedom in 
SFR calibration, (b) Freedom of cr^ = 0.2 in SFR calibration (see 
explanation in text), (c) Freedom of cr^ = 0.5 in SFR calibration. 



The horizontal lines represent the three values of the back- 
ground derived in § 13.11 The rounded, vertically extended 
contours represent joint confidence intervals for these two 
parameters. We construct these by multiplying the probabil- 
ity distributions of the two parameters together, and finding 
the contours that enclose 65%, 95%, and 99.8% of the prob- 
ability. The "squashed" appearance of the contours comes 
from the non-Gaussian probability distribution of the back- 
ground light (see H3.1l for an explanation of how we estimate 
its distribution). 

The diagonal contours in panel (a) depict the 65%, 95%, 
and 99% 2-dimensional confidence regions derived from the 
ensemble of star formation histories. These contours form a 
narrow band. If we consider only one dimension, our com- 
puted values of Jebl are in reasonable agreement with ob- 
servations, while the stellar mass is rather high as mentioned 
before. The main point, however, is that the two sets of con- 
tours barely overlap; one needs to go out to the Sa contours 
to find agr eement. We obtain essentially identical results if 
we use the lCole et alj i2001f) parametric form for the star 
formation history instead of the form in Equation |S] 

The remaining two panels show the results for samples 
with increasing levels of assumed systematic error in the 
overall calibration of the SFRs, i.e. a random dispersion in 
the calibration of p-kform- This greatly expands the allowed 
joint confidence region of AT-band hght and EBL so that, 
as we just discussed, the former is no longer inconsistent 
with the observations. However, the confidence region still 
forms an extremely narrow band, which overlaps the obser- 
vational contours at only the 2 or 3(t level. Even with this 
increased freedom, one can find SFR histories that agree well 
with either the observed K light or the observed bolometric 
background, but not with both quantities simultaneously. 

We consider other proxies for the stellar mass as well. 
Plotting Jebl versus pj (instead of pk) in Figure |S| we see 
that the 2cr contours still do not overlap. If instead we use 
the SDSS z-band density, p^, the contours still do not sig- 
nificantly overlap. The slightly better agreement in 2:-band 
might suggest that there may indeed be a small bias in 
the 2MASS surveys owing to missing low surface brightness 
galaxies. However, our extinction estimates become increas- 
ingly suspect as we go to these shorter wavelengths. Finally, 
we plot the contours of stellar mass itself. With the larger ob- 
servational uncertainties in this quantity, the disagreement is 
now at less than the 2a level, though the impression remains 
that either the observed background is too high or the ob- 
served mass too low. We regard the K-ha,nA plots as more 
fundamental, since they plot a directly observed quantity 
and one that is tied as closely as possible to the dominant 
stellar mass range. 

In essence, we find that the bolometric background is 
too bright relative to the present day _ft'-band hght density, 
given a Salpeter-like IMF and the shape of the SFH implied 
by our fits to the observations. Changes to the normalisation 
of the SFH move predictions along the diagonal in Figures |K| 
and|^and, therefore, do nothing to resolve this discrepancy. 
Our (7A = 0.5 models allow substantial variations in the SFH 
shape as well as normalisation, but these have limited power 
to solve the problem. The ratio between the background pro- 
duction function Sebl{z) and the i^-band luminosity Pk{z) 
has a maximum at z « 1 and falls off slowly at higher and 
lower redshifts. Hence, by the mean value theorem, the ra- 
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Figure 6. The contours of EBL vs. (a) J-band stellar light 
density, (b) 2-band stellar light density, and (c) stellar mass den- 
sity allowed by the range of SFHs. Except for substituting these 
quantities for i^T-band light, the plots are the same as in the previ- 
ous figure. In panel (c), the asterisk shows the estimated relation 
of stellar mass and EBL from Madau ct al. (2001), rcscalcd to 
our IMF, and the cross shows the result from the simulation of 
[Nagaminc et al. |200C|) (see text for details). 



tio between the EBL and the emitted A'-band luminosity 
density cannot be larger than 



/ Jebl{z 



V Piiiz) 



_ I Sebl{z 
= 1.21 X 10"^bgu[LQ^ Mpc 



max ' """ 

4-Kf '" 



for a Salpeter or diet Salpeter IMF, including the extinc- 
tion correction. We plot the maximum ratio as the diagonal 
dotted line in Figure |K| Because much of the star forma- 
tion in our calculations takes place at 2; ~ 1, the predicted 
ratio for our best-fit SFH, 9.0 x lO"*bgu[L0^ Mpc"^]~\ 
is already 75% of this upper limit, so even a 5- function 
SFH a.t z — 1 can produce only a modest increase. The 
ratio of our middle trace EBL to our best estimate pK 



1.8 X 10"^bgu[LQ^ Mpc" 



Getting down to 1.0 



10~'^bgu[LQjy Mpc""^]"^ requires taking our minimal esti- 
mate of 50 bgu for Jebl (the lower trace in Figure and 
increasing our averaged value of pK ~ 4.2 x 10* Lq^ Mpc~'^ 
by 15%. 

With our best SFH fit, 50% of the background is pro- 
duced in stars that form at z < 1.0, and the central 
80% comes from the stars that form in the redshift range 
0.27 < z < 2.6. This can be compared with the present-day 
stellar mass, for which 50% comes from z < 1.2, and the 
central 80% comes from the redshift range 0.33 < z < 3.2. 
For the A'-band luminosity density, 50% comes from stars 
formed at z < 0.9, and the central 80% comes from the red- 
shift range 0.1 < z < 2.7. Therefore, the EBL principally 
comes from stars formed at the same redshifts as those that 
dominate the current mass or K-hand fight, which explains 
why the correlation is so tight. The background energy is 
emitted at still lower redshifts. Again, using our best esti- 
mate of the SFH, 50% comes from z < 0.65, and the central 
80% comes from the redshift range 0.11 < z < 1.9. 

Even if we consider changes in the star formation 
history that are more radical than those allowed by our 
random-calibration formalism, it is still difficult to change 
the ratio of background light to stellar mass significantly. 
The "fossil re cord" SFH, estimated from the spectra of lo- 
cal galaxies bv lHeavens et al.N2004) . peaks at a significantly 
lower redshift than our typical fits, which is more favourable 
to the creation of EBL. However, the ratio of A'-band light 
to stellar mass, using this star formation history, would rise 
by only a few percent and the ratio of EBL to A'-band lu- 
minosity would actually drop by a few percent, since stars 
at lower redshifts are even more capable of producing K- 
band light. Two independent calculations of the EBL versus 
stellar mass from the literature, which we discuss below, are 
also shown in Figure |S| These values bolster our conclusion 
that the ratio is insensitive to reasonable changes to the 
SFH. Since all allowable star formation histories imply that 
most stars are old (Figure]^, even significant shifts in the 
distribution of their formation redshifts only results in mi- 
nor changes to the majority of the stellar ages and to the 
resulting stellar energy output. 

The values of p* and pK calculated from the observa- 
tions are proportional to the value of the Hubble constant. 
Ho- However, the observed EBL and pi, form, integrated from 
the SFH, are both independent of the Hubble constant, as 
is F-t, also computed from the integrated ptform owing to 
the invariance of pbaryon derived from microwave background 
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measurements. So an uncertainty of 7kms~^ Mpc~^ in our 
canonical value of Hq = 70kms~^ Mpc~^ will add a 10% 
uncertainty when comparing p* to these other quantities. 
(The changes from offsets in the lookback-time arguments 
in Equations 1-3 are smaller still and can be neglected.) 

In Figure |S1 we have not taken into account the sys- 
tematic error from uncertainties in the Hubble constant or 
other systematics such as uncertainties in the stellar popula- 
tion tracks and synthesis, the stellar metallicity, or the mean 
NIR extinction. Taken together, these systematics may com- 
bine to a net 25% error, but we do not place much faith in 
this estimate. With systematics in mind, an observational 
solution to the background problem seems reasonable, per- 
haps giving values close to pK = 5.0 x 10^ Mq Mpc~^ and 
Jebl = 55 bgu. We emphasise that for the background light 
to be so near its lower limit, a large set of absolute back- 
ground detections would need to be wrong, and if future 
observations substantiate the higher detections in any re- 
gion of the spectrum, an "observational" solution will be 
essentially ruled out. 

O ur conclusion of a p o tentia l con flict runs contrary t o 
that of'Madau fc Pozzettil i2000l) and lMadau et al.l J200ll) . 
who find very good agreement between the background ra- 
diation and the stellar density for their assumed star for- 
mation histories. '^ These authors use a stellar density of 
p* = 6.8 X 10®/i Mq Mpc~^ based on B -band surveys, and 
find a best EBL value of Jebl = 55 bgu llMadau fc Pozzettil 
I2OOOI) or 60 bgu (Madau et al.ll200ll) . We plot this result in 
Figure |S| after translating to our IMF and including mass 
recycling, and it is in good agreement with our ratio of Jebl 
to pi,. Since that time, the estimates of p^hT^ , using larger 
and more accurate NIR-based galaxy surveys, have dropped 
by nearly a factor two, a change only partially mitigated by 
the change in the value of the Hubble constant from their 
assumed value of 50 kms~^ Mpc~^ to our assumed value 
of 70 kms~^ Mpc~^. Also, additional observations of the 
background light have raised the most likely value of Jebl, 
making the previous "best values" closer to our lower limit. 
These two changes explain their discrepancy with our con- 
clusions.^ 

■Nagamine et al.l (|2000J) computed the star formation 
history using an LCDM hydrodynamic simulation, and find 
from this history that Jebl ~ 81 bgu and that the stellar 
density, using a Salpeter IMF, is f2* — 0.012. Translating 
this to our diet Salpeter IMF and including mass recycling 
implies p* — 8.0 x 10*MqMpc~^, over twice as large as 
what is acceptable today. We also show this result in Fig- 
ure |n] Nagamine et al.'s ratio of Jebl to p* is only 15% 
lower than ours, and this slight difference probably owes 
to their high level of star formation at late times, which 
exceeds observational constraints. In a more recent paper 
(which appeared in preprint form as we were finali sing this 
manuscript for submission) , iNagamine et al.l (|2D06|) present 
empirical and numerical models of the cosmic SFH which. 



^ The int egrated background from the "realistic" star formation 
history in lMadau fc Pozzettil 120001) appears to be incorrect, but 
this i s corrected in the later conference proceedings jMadau et alj 
|2001^ along with slight changes to the assumed star formation 
history. The ratio of background radiation to stellar mass density 
for the cases considered in the latter paper is then quite insensitive 
to the star formation history, in agreement with our results above. 



they argue, are consistent with both the extragalactic back- 
ground light and the local bolometric luminosity density. 
However, they adopt a background value of Jebl ~ 45 bgu, 
which is even lower than the value of 50 bgu implied by our 
lowest trace through the data in Figure^ in essence, they as- 
sume that the integrated count estimates of the background 
are correct and that all of the absolute measurements are 
incorrect. We also regard the local 7^-band light density as 
a more robust "fossil" constraint on the SFH than the local 
bolometric light density, because the latter is more difficult 
to estimate observationally and is more sensitive to contri- 
butions f rom very rec ent star formation. We concur with 
Nagamin e et al.l i2006r that the directly estimated SFH, ex- 
tragalactic background light, and local light density can be 
reconciled if one takes a low value for the second and a high 
value for the third, but this solution requires pushing the 
systematic uncertainties of the observations to their limits. 



5 SOME POSSIBLE SOLUTIONS TO THE 
BACKGROUND PROBLEM 

With our assumed diet Salpeter IMF, the local stellar mass 
density (or light density) appears to be inconsistent with 
the directly estimated cosmic star formation history and 
the extragalactic background light. The local light suggests 
a stellar baryon fraction Kt ~ 0.06, while the integrated 
star formation history suggests a higher value of F-^ « 0.09. 
The EBL suggests values that are similar or even larger. 
Ft > 0.08 — 0.20, and the fairly hard lower limit on the ob- 
served EBL makes it difficult to accommodate without mak- 
ing the stellar mass density too high. The calibration of the 
SFH is uncertain enough that it can be made consistent 
with either of the other two constraints individually, but 
not simultaneously. Systematic uncertainties in the absolute 
background measurements and the completeness of the lo- 
cal galaxy census leave room for an "observational" solution 
to this conflict. However, if future measurements accurately 
constrain the stellar mass density and background light to 
the most probable values we have obtained here, there will 
be a substantial discrepancy. 

Let us, therefore, consider the alternative solution of 
extra energy in the background beyond that provided by 
a Salpeter IMF. Solving the background problem using an 
additional source of energy requires an extra 10-40 bgu, or 
0.4-1.7 X lO"^'' ergcm"^, which is comparable to the energy 
density of stellar light from all the detected galaxies, even 
discounting the effect of redshift. This amounts to ~ 10~^ of 
the closure energy density, or about 12-50 keV per baryon. 
Thi s is a steep requirement f or any possible energy source 
fsee lFukugita fc Peeblea2004l for a useful list of candidates) . 
For exotic sources, the requirement that the excess appear 
in the optical/NIR or FIR portions of the spectrum is an 
additional hurdle. 

The energy requirement alone immediately rules out 
most possibilities, including supernova-driven shocks and 
gravitationally powered cooling radiation. Decaying cosmic 
neutrinos or decaying dark matter would have sufficient en- 
ergy, but producing the required decay rate and appropriate 
photon energies would require substantial fine-tuning. 

AGN are a known contributor to the background at the 
~ 2 bgu level, as discussed in § 13.11 A large (e.g., factor 
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of 5-10) increase in this contribution is possible only if the 
local census of fossil black holes is seriously incomplete, per- 
haps because the great majority of supermassive black holes 
have been ejected from galaxy centres in mergers. A num- 
ber of authors have suggested Pop III stars or their black 
hole remnants as sources of the high background measured 
in the l — 4jj, region, but lMadau fc Silkl i2005f) show that the 
energetic and chemical constraints on such a population are 
difficult to satisfy; furthermore, they have been suggested 
to explain a set of measurements that we are already not 
including in our fiducial estimate of Jebl- 

Thus, non-stellar solutions to the background problem 
do not appear promising. On the other hand, we have not 
yet considered modifications of the IMF. An IMF that is 
biased to intermediate- or high-mass stars, compared to the 
diet Salpeter IMF we have used up until now, could have a 
significant effect on the stellar emission. These stars would 
mostly burn out by the present day, thus contributing to 
the observed SFH and EBL but not to the A'-band light 
density. As noted in the introduction, improving agreement 
among these observations requires changing the IMF above 
the turnoff mass 0.9Mq of old stellar populations; changes 
below this mass renormalise all three quantities by the same 
factor.'^ 

In Figure [7^ we show several IMF s that have been pro- 
posed in the literature ( Salpeter 195 51: iMiller fc Scaldll979l: 
Kcnnicutt 1983; Kroupa ct al. 1993; Krouoa 2001; Chabrierl 
2003 : Baldrv fc Glazebrook 2003) . all normalised to the 
same value at 1 Mq . We also show th e approximate a llowed 
region of the IMF as determined by iKroupal (|2002) ■ This 
allowed region was established from studies using several 
different techniques covering different mass ranges; hence 
the normalisations above and below 1 Mq are poorly con- 
strained relative to one another. The slopes from 1 to 15 Mq 
in the plot are 1.7 for Kroupa (1993), 1.6 for Miller-Scalo 
(1979), 1.5 for Kennicutt (1983), 1.35 for Salpeter (1955), 
1.3 for Kroupa (2001, mean Galactic- field form in Equation 
2) and Chabrier (2003), and 1.2 for iBaldrv fc Glazebrook 
J2003h IMFs. 

Now let us repeat the exercise leading to Figure]^ but 
using the different IMFs. We normalise the star formation 
histories by the amount of UV light produced in a 10* yr 
burst. (Using the bolometric luminosity at the same burst 
length would give a similar scaling.) Figure |7|d shows the 
resulting K-hanA-JEBL contours. There is a single overall 
trend in this plot. The stars that dominate the ii'-band 
light today (pk), the EBL, and the UV light are of pro- 
gressively higher mass, so as the IMF slope above 1 Mq 
is made shallower, the contours shift to lower values of pK 
and Jebl, while at the same time the ratio of Jebl/ pK in- 
creases. (Changes to the IMF below 1 Mq have little effect 
on the A'-band or bolometric light.) Because pK is most 
likely better known than the calibration of the SFH, mak- 
ing the IMF shallower above 1 Mq actually raises the most 
plausible background value, despite the resulting decline in 
the mean Jebl of the contours. 



^ iNagamine et all |200g) conclude that IMF changes have little 
impact on the predicted EBL, but that is because they only com- 
pare two cases (Salpeter and Chabrier) that arc nearly identical 
above 1 Mq . 



We consider the Salpeter IMF to be the point at which 
the agreement is marginally acceptable, and rule out IMF 
slopes steeper than the Salpeter value of —1.35. This ex- 
cludes a n umber of IMFs c ommonly used in the literature. 
The Baldrv fc Glazebrookl JJOOS) IMF is the most accept- 
able; these authors estimate a high-mass slope of — 1.15±0.2, 
which independently rules out the steepest IMFs. Alone 
among the IMFs in Figure [TJa, the Baldry-Glazebrook IMF 
was derived from the global SFH, in this case by comparing 
the spectrum of the local luminosity density to observation- 
ally plausible shapes of the SFH. It is interesting that it 
agrees best with our independent test, perhaps indicating 
that the estimates from Galactic data are somehow biased 
or that the Galactic IMF differs from the Universal average. 

However, none of the standard IMFs from the literature 
give more than marginal agreement with our joint contours 
in pK and Jebl- This will pose a problem if future obser- 
vations constrain the background to a high value. We now 
consider several examples to see what is required to get sig- 
nificantly higher levels of background light. 

We first consider the case of universal IMFs. For our 
purposes this does not necessarily imply a completely in- 
variant IMF from one galaxy to another, but does require 
that there be no systematic change with redshift. It is use- 
ful to distinguish between a top-heavy IMF, rich in high- 
mass (M > 8Mq) stars relative to Sun- like stars, and a 
middle-heavy or "paunchy" IMF, rich in intermediate-mass 
(IMq < M < 8Mq) stars. As discussed above, the par- 
ticular range that would be most efficient in boosting Jebl 
relative to the UV/FIR or K-h&nA light is 1.5-4 Mq. 

In Figure [TJ; and Tabled we introduce three additional 
IMFs. The first is an "extreme top-heavy" IMF (the slope 
here is 0.95, so that the mass converges only at the upper 
end of the mass function). As an example of a "paunchy" 
IMF, we place a break at 4Mq, and for a more extreme or 
"obese" example, we place a discontinuity at 1.5 Mq. For 
comparison, we also include a diet Salpeter IMF, and we 
use the Baldry-Glazebrook IMF as an example of a modestly 
top-heavy IMF. 

We plot results in the pk-Jebl plane for these IMFs in 
Figure |7|i. Here, we only represent the results in a schematic 
way to avoid overlapping contours. The top-heavy examples 
continue the general trend down and to the left in the plot. 
The paunchy examples, by contrast, go only to the left since 
they have enhanced bolometric emission relative to the K- 
band or UV light. 

Our paunchy and extreme top-heavy IMFs resid e within 
the hatched allowed regions of lKroupal (|2001() and iKroupal 
(2002). For 0.5-1.0 Mq, the paunchy IMF is in marginal 
conflict with Kroupa's limits on the slope, but this is based 
on only a few samples with large error bars and potentially 
large systematics. The Baldry-Glazebrook and extreme top- 
heavy IMFs are in good agreement with all of these limits. 
The obese IMF, however, is too "paunchy" to be consistent 
with Milky Way observations. 

The need to join together estimates in different mass 
ranges from different environments into a single IMF con- 
tributes to an uncertainty in the observed shape of the IMF 
dMasseviiiggS : Kroupa 2002 ; Chabrier 2003: Lada fc Ladal 
l2003ll . As shown in Figure 5 of lKroupal ^l2002^ ■ estimates of 
the slope in the range 0.8-3 Mq show a particularly large 
scatter, because neither the Galactic field star or young clus- 



© 0000 RAS, MNRAS 000, 000-000 



18 Mark A. Fardal et al. 



10.00 



1.00 



0.10 



0.01 L 
0.1 




Kroupa 1993 
Miller-Scalo 

■ Keniiicutt 

■ Diet Salpeter 

■ Kroupa 2001 

■ Chabrier 

■ Baldry-Glazebrook 



1.0 10.0 

Mass (Mgl 



100.0 




K-band light density (L^j Mpc ) 



10.00 



1.00- 



0.10 



0.01 
0.1 




■ Diet Salpeter 

■ Baldry-Glazebrook 

■ Extreme top-heavy 

■ Paunchy 

■ Obese 

■ Avg. of Variable 




1.0 10.0 

Mass (Mg) 



100.0 



20 



10' 



ADiet Salpeter 
^Baldry-Glazebrook 
Paunchy 

Obese 

▼Extreme top-heavy 



rVariable 



10' 



K-band light density (Lg; Mpc ') 



Figure 7. (a) The standard IMFs considered in the paper. Th e vertical axis i s the mass contained in the IMF per unit logarithmic 
stellar mass. The shaded region is the allowed region taken from iKroupal i2002^ . (b) The same as Figure 13 but showing the effects of 
using the different IMFs in the previous panel. We scale each SFH to account for the UV luminosity in a 10* yr burst for that particular 
IMF. The contours reflect a systematic error of a a = 0.2 in the SFH calibration, (c) Shows the new IMFs introduced in the paper. We 
also plot Baldry-Glazebrook as an example of a moderately top-heavy IMF and Diet Salpeter for comparison, (d) Like panel (b) but for 
the IMFs shown in panel (c). To avoid the confusion of overlapping contours, we only show the values from the best fit and indicate the 
approximate contour extent by a dashed line. 



ter methods is very effective in this region. Kroupa notes 
that the low and high mass regions of the IMF are observa- 
tionally almost disjoint, owing to the lack of overlap between 
the mass ranges of the most accurate methods. Thus if any- 
thing, the hatched regions in Figures |7Ji and|7J; are too re- 
strictive in the region aroun d 1 Mr:^ , m aking a paunchy IMF 
more plausible. For example. ISirianni e t al. (2000) used HST 
photometry of the large star cluster 30 Dor in the LMC to 
obtain isochrones in the colour-magnitude diagram. Their 
resulting IMF has a slope of 1.28 ± 0.03 for A/ > 2.1 Mq, 
and 0.27 ± 0.08 for M < 2.1 Mp-., at least down to 1.3 Mp 



^0 



where the noise becomes large. This is tentative support for 
a paunchy IMF, at least in some active star-forming envi- 
ronments, which may be more representative of the main 
contributors to the cosmic SFH. 



Interestingly, iKroupal l|200J) finds evidence for a larger 
number of low-mass stars in the present-day IMF than in 
the mean Galactic IMF, which is dominated by stars formed 
within the last 5 Gyr. This is tentative evidence that the IMF 
has been increasingly weighted to stars above 1 Mq in the 
past, and it would not be surprising if the trend continued to 
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the global population of stars, which is still older on average 
than the local Galactic disk. 

Let us now consider an IMF that varies systematically 
with cosmic time. There are both theoretical and obser- 
vational motivations for considering a top-heavy IMF for 
high-redshift or high-SFR galaxies. Problems that this has 
been invoked to solve include gas consumption timescales 
in starbursts (Rickc, Lokcn, Ricko, & Tamblvnl ll993h . sys- 
tema tic variations of m etallicity and M/L ratios in ellipti- 
cals lZeBf_fc_SillJL99a), and large iron abundan ces in clus- 
ters l|zepf fc Silklll996l : iMaoz fc Gal-Yamll2004) . It has of- 
ten been suggested that in extreme star formation envi- 
ronments, the formation of low-mass stars should be sup- 
pressed either by the global increased temperature of the 
medium, which basically rais es the Jeans mass, or by feed - 
back from high-mass stars JMassevI Il998l : iLarsonl Il99d) . 
However, proving this by observing low and high mass stars 
together in the same environment is very difficult. Simu- 
lations suggest that primordial metal- free stars, at least, 
form with an ex tremely top-heavy IMF (with a minimum 
mass ~ 100 Mg: lAbel. Brvan. fc Normanli2002h . while sim- 
ulations of starburst ga laxies also suggest a top-heavy IMF 
JPadoan. Nordlund. fc Jones 1997) . The Milky Way IMF 
shows two cha racteristic mass scales, at ~0.01 and ~0.5 Mq 
^Kroupal2001^ . whose origin is not well understood. It would 
not be surprising if these characteristic masses shifted with 
redshift, metallicity, or SFR. 

As a toy model to investigate this idea, we suppose 
that there are both "normal" and "rapid" modes of star 
formation and associate the latter with starbursting galax- 
ies, rnaking it much more important at high redshifts (cf. 
iLaeache et al J 120031) . We assume both modes have a Ken- 
nicutt IMF, but sharply truncate it in the "rapid" mode 
below 3Mq. We assume that the normal mode contributes 
a fixed amount of star formation at all redshifts, which is 
90% of the star formation rate today, and that the rapid 
mode contributes the remainder. In this model the rapid 
mode quickly becomes the dominant form as the redshift 
increases. At very high redshifts, where the SFR rate falls 
below 90% of the rate today, we just assume that all the 
star formation is in the normal mode. Each component is 
scaled according to its UV luminosity at 10* yr to match 
the specified fraction of the best-fit SFH. 

We show the time-averaged IMF produced by this toy 
model in FigureEt, and its resulting location in the pk-Jebl 
plane in Figure |7|l. Clearly, this type of model can produce 
large background levels, especially if the calibration of the 
SFH is poorly known at high redshift. We note that from 
an energetic standpoint this is an inefficient solution, since 
much of the large energy output at early times is redshifted 
away. It would be more efficient to make the top-heavy mode 
dominant at low redshift, but this would then lead to con- 
flicts with local observational constraints on the IMF. This 
particular model can be described as paunchy, since it has 
a sudden jump at 3 Mq . 

We conclude that a paunchy or top-heavy IMF is ob- 
servationally allowed, and might indeed be required to solve 
the background problem, depending on how future observa- 
tions sort themselves out. The consequences of a top-heavy 
IMF have been frequently considered. A paunchy IMF also 
has numerous consequences, which are similar in some ways 
and different in others to a top-heavy IMF. For a given SFH, 



it increases the number of F, A, and late B main-sequence 
stars and the number of bright giants. This will increase 
the optical brightness of galaxies relative to the IR, or raise 
the specific SFR that is inferred assuming a standard IMF. 
Currently, the SFH deduced from the fossil record is signif- 
icantly higher than direct SFH measurements at lookback 
times of ^ 5 Gyr, and signifi cantly lower at < 1 Gyr and 
>8Gyr JHeavens et alj|2004ln . This suggests that there are 
fewer low-mass and high-mass stars than one expects from 
simply integrating the observed SFH, assuming a standard 
IMF. Assuming a paunchy IMF would help solve this dis- 
crepancy. 

With a paunchy IMF, post-starburst "E-f A" galaxies 
will appear more frequently. The cosmic luminosity density 
will also be mo dified, and it would be worthwhile to repeat 
the analysis of iBaldrv fc Glazebrookl i2003ll . with a break 
at intermediate stellar masses, to see whether the high-mass 
slope we have assumed for the paunchy form could be con- 
strained. The passive evolution of galaxies would appear 
more rapid with a paunchy IMF, so that the apparent mass 
function of high-redshift red-sequence galaxies would shift 
to higher masses than if one assumed a standard IMF. The 
intermediate-mass stars will leave large numbers of white 
dwarfs, amounting to ~ 20% of the total stellar mass, which 
could be constrained by microlensing surveys. Finally, the 
chemical byproducts of intermediate-mass stars are very 
different from those of high-mass stars, tending to make 
"a-suppressed" abundances. Hence, the analysis of cosmic 
abundances may be able to place limits on the shape of this 
type of IMF. 



6 CONCLUSIONS 

We have investigated the relationship between the total 
amount of extragalactic background light, the cosmic star 
formation history, and the stellar population today. We find 
that there is tension, and arguably outright conflict, between 
the background light and NIR surveys that probe the local 
stellar mass density around 1 Mq , if we assume a standard 
Salpeter IMF slope. We have paid particular attention to 
uncertainty in the star formation history; we find that it 
contributes only a small dispersion to the predicted ratio 
of background light and NIR light, even though the total 
amount of star formation is highly uncertain. 

This ratio, however, is IMF-depende nt. Universal IMF s 
with steeper slopes a bov e IM^, such as | Kennicutt] il983l) , 
iMiUer fc Scald Jl97gl) . or lKroupa et alj ||l99^~areniled out 
by the current minimum energy requirements. If the higher 
absolute background measurements should be substantiated 
in the future, it conversely implies that the average IMF is 
slightly top-heavy, or more likely (given limits on the ob- 
served star formation history), middle-heavy or "paunchy", 
i.e. rich in 1.5-4 Mq stars compared to the Salpeter IMF. 
Alternative solutions in which non-stellar sources make large 
contributions to the background light appear unlikely. 

We note that every quantity that one can measure for 
the stellar population of galaxies is IMF-dependent, and 
many secondary conclusions drawn from the study of galac- 
tic stellar populations, e.g. galaxy merger rates or rates of 
gas infall, depend in an indirect way on the IMF. The ex- 
treme difficulty of deriving the stellar IMF from first prin- 
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ciples, the fundamental difBculty of establishing the IMF 
above 1 Mq from local observations, and the theoretical sus- 
picion that the IMF varies with environment all suggest that 
galactic astronomers may have to deduce the form of the 
IMF from their own data. For that, we will need to develop 
good tests involving closed systems of gas and stars. Since 
the Universe is the ultimate closed system, it seems likely 
that studies of cosmic- averaged quantities will be crucial to 
establish the form of the IMF. The study of cosmic star for- 
mation and background light in this paper is a step in that 
direction. 
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Panter, and Steve Schneider for useful conversations. This 
research was supported by NASA grants NAG5-13102, 
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